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Abstract — Additive codes over GF(9) that are self-dual with 
respect to the Hermitian trace inner product have a natural ap- 
plication in quantum information theory, where they correspond 
to ternary quantum error-correcting codes. However, these codes 
have so far received far less interest from coding theorists than 
self -dual additive codes over GF(4), which correspond to binary 
quantum codes. Self-dual additive codes over GF(9) have been 
classified up to length 8, and in this paper we extend the complete 
classification to codes of length 9 and 10. The classification is 
obtained by using a new algorithm that combines two graph 
representations of self-dual additive codes. The search space is 
first reduced by the fact that every code can be mapped to a 
weighted graph, and a different graph is then introduced that 
transforms the problem of code equivalence into a problem of 
graph isomorphism. By an extension technique, we are able to 
classify all optimal codes of length 11 and 12. There are 56 005 876 
(11, 3 11 , 5) codes and 6493 (12, 3 12 , 6) codes. We also find the 
smallest codes with trivial automorphism group. 

Index Terms — Self-dual codes, additive codes, codes over 
GF(9), graph theory, classification, nonbinary quantum codes. 

I. Introduction 

ADDITIVE codes over GF(9) of length n are GF(3)- 
linear subgroups of GF(9) n . Such an additive code 
contains 3 fe codewords for some < k < 2n, and is 
called an (n, 3 fe ) code. A code C can be defined by a k x n 
generator matrix with entries from GF(9) whose rows span C 
additively. We denote GF(9) = {0, l,w,w 2 , . . . ,uj 7 }, where 
uj 2 =u + l. Conjugation of x 6 GF(9) is defined by 
x = x 3 . The trace map, Tr : GF(9) -)• GF(3), is defined 
by Tr(x) = x + x. Following Nebe, Rains, and Sloane (TJ, 
we define the Hermitian trace inner product of two vectors 
u,v e GF(9)" by 

(it, v) = u> 2 (u ■ v — u ■ v) = Tr(w 2 M • v), 

where multiplication by uj 2 is necessary because the skew- 
symmetric bilinear form (u-v — u-v) does not take values in 
GF(3) [ 1 1. We define the dual of the code C with respect to the 
Hermitian trace inner product, C = {u e GF(9) n | (it, c) 
f) for all ceC}. C is self-orthogonal if C C C x . If C = C^, 
then C is self-dual and must be an (n, 3") code. The class 
of trace-Hermitian self-dual additive codes over GF(9) is also 
known as 9 H+ (TJ. The Hamming weight of u, denoted wt(tt), 
is the number of non-zero components of it. The Hamming 
distance between u and v is wt(u—v). The minimum distance 
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of the code C is the minimal Hamming distance between any 
two distinct codewords of C. Since C is an additive code, 
the minimum distance is also given by the smallest non-zero 
weight of any codeword in C. A code with minimum distance d 
is called an (n, 3 fc , d) code. The weight distribution of the code 
C is the sequence (Aq, A\, . . . , A n ), where A4 is the number 
of codewords of weight i. The weight enumerator of C is the 
polynomial 

n 

w{x, y ) = Y. A * xn ~ l y l i 

where we will denote W{y) — W(l,y). It follows from 
the Singleton bound that any self-dual additive code 
must satisfy d < [§J + 1. A code is called extremal if it 
has minimum distance |_§J + 1, and near-extremal if it has 
minimum distance [§J- If a code has the highest possible 
minimum distance for the given length, it is called optimal. A 
tighter bound exists for codes over GF(4) J3], but in general 
the Singleton bound is the best known upper bound. Codes that 
satisfy the Singleton bound with equality are also known as 
maximum distance separable (MDS) codes. The well-known 
MDS conjecture implies that self-dual additive MDS codes 
over GF(9) must have length n < 10. We have shown in 
previous work |4| that there are only three non-trivial MDS 
codes, with parameters (4,3 4 ,3), (6,3 6 ,4), and (10, 3 10 , 6), 
given that the MDS conjecture holds. 

Two self-dual additive codes over GF(9) are equivalent if 
the codewords of one can be mapped onto the codewords of 
the other by a transformation that preserves the properties of 
the code, i.e., weight enumerator, additivity, and self-duality. 
It was shown by Rains [2 | that this group of transformations 
is Sp 2 (3) I Sym(n), i.e., permutations of the coordinates 
combined with operations from the symplectic group Sp 2 (3) 
applied independently to each coordinate. Global conjugation 
of all coordinates will also preserve the properties of the code, 
and codes related by this operation are called weakly equiv- 
alent In this paper, we classify codes up to equivalence, 
i.e., we do not consider global conjugation. Let an element 
a + buj e GF(9), be represented as (£) e GF(3) 2 . We can 
then premultiply this element by a 2 x 2 matrix. The group 

**>-<(! 0) 

has order 24 and contains all 2 x 2 matrices with elements from 
GF(3) and determinant one. The order of Sp 2 (3) I Sym(n) is 
24™n!, and hence this is the total number of maps that take a 
self-dual additive code over GF(9) to an equivalent code j2). 
By translating the action of Sp 2 (3) on (?) into operations on 
elements c = a + buj G GF(9), we find that the operations 
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we can apply to all elements in a coordinate of a code are 
c M> xc if x 4 = 1, or c 1-4 xc if x 4 = — 1, given x G GF(9), 
and a + bu i->- a + yb + bu, given y G GF(3). 

A transformation that maps C to itself is called an au- 
tomorphism of C. All automorphisms of C make up the 
automorphism group of C, denoted Aut(C). The number of 
distinct codes equivalent to C is then given by rs^t^TI • ^ ne 
equivalence class of C contains all distinct codes that are 
equivalent to C. By adding the sizes of all equivalence classes 
of codes of length n, we find the total number of distinct codes 
of length n, denoted T n . The number T n is also given by a 
mass formula which was described by Hohn [5] for self-dual 
additive codes over GF(4) and is easily generalized to GF(9): 



T n = n(3 i 

1=1 



24™n! 



T[ |Aut(C,; 



(1) 



where t n is the number of equivalence classes of codes of 
length n, and Cj is a representative from each equivalence 
class. The smallest possible automorphism group, called the 
trivial automorphism group, of a self-dual additive code over 
GF(9) is {/, —I}, i.e., it consists of global multiplication of 
coordinates by 1 or — 1. By assuming that all codes of length 
n have a trivial automorphism group, we obtain from the mass 
formula a lower bound on t n , the total number of inequivalent 
codes. 



tn > 



24"n! 



(2) 



Note that when n is large, most codes have a trivial auto- 
morphism group, so the tightness of the bound increases with 
n. As we will see in Section I VIII for n = 10, 80% of all 
codes have a trivial automorphism group, and the bound (01 
underestimates t%o by just 19%. 

Any linear code over GF(9) that is self-dual with respect 
to the Hermitian inner product, (u, v) = u ■ v, is also a self- 
dual additive code with respect to the Hermitian trace inner 
product. The class of Hermitian self-dual linear codes over 
GF(9) is also known as 9 H (T). The operations that map a 
self-dual linear code to an equivalent code are more restrictive 
than for additive codes, since GF(9)-linearity must now be 
preserved. Only coordinate permutations and multiplication of 
single coordinates by x G GF(9) where x 4 = 1 is allowed. It 
follows that only additive codes that satisfy certain constraints 
can be equivalent to linear codes. Such constraints for codes 
over GF(4) were described by Van den Nest |6] and by Glynn 
et al. 0. An obvious constraint is that all coefficients of 
the weight enumerator, except Aq, of a linear code must be 
divisible by 8, whereas for an additive code they need only be 
divisible by 2. To our knowledge, no complete classification 
of Hermitian self-dual linear codes over GF(9) have appeared 
so far, but several authors have studied this class of codes 
and suggested a number of constructions [8|-[11|. Checking 
whether a self-dual additive code over GF(9) is equivalent 
to a linear code is non-trivial, since there are 6™ coordinate 
transformations in Sp 2 (3)™ that could transform a non-linear 
code into a linear code. Our classification of self-dual additive 
codes could be a useful starting point for also studying linear 
codes, but this is left as a problem for future work. 



Trace-Hermitian self-dual additive codes over GF(q) exist 
for q = m 2 , where m is a prime power JT|, and the class of 
self-dual additive codes over GF(q) is called q H+ . The first 
case, 4 H+ , has been studied in detail, in particular since an 
application to quantum error-correction was discovered ||3]. 
We have previously classified self-dual additive codes over 
GF(4) up to length 12 JT2). Self-dual linear codes over GF(4) 
have been classified up to length 16 lfT3l by Conway, Pless, 
and Sloane. This classification was recently extended to length 
18 IH and 20 US) by Harada et al. The next class of self-dual 
additive codes, 9 H+ , has received less attention, although these 
codes have similar application in quantum error-correction 0, 
lfl6ll . where they correspond to ternary quantum codes. We 
have previously classified self-dual additive codes over GF(9) 
up to length 8 [4], as well as self-dual additive codes over 
GF(16) and GF(25) up to length 6. Another type of self-dual 
code over GF(9) is known as 9 E (H and is self-dual with 
respect to the Euclidean inner product, (u, v) = u ■ v. There 
is no additive variant of these codes, and this family will not be 
considered in this paper. Again, some constructions have been 
described |10|, but no complete classifications of Euclidean 
self-dual codes over GF(9) have been given. 

In Section [n] we briefly review the connection between 
trace-Hermitian self-dual additive codes and weighted graphs. 
An algorithm for checking equivalence of self-dual additive 
codes over GF(9), which is a generalization of a known 
algorithm for linear codes ifTTIl . is described in Sections [Til] 
and [IV] Combining this algorithm with the weighted graph 
representation, and some other optimizations, enables us to 
classify all self-dual additive codes over GF(9) of length up 
to 10 in Section [V] In particular, all near-extremal codes of 
length 9 and 10 are classified for the first time. We also find 
the smallest codes with trivial automorphism group. Using an 
extension technique described in Section [VTl we are then able 
to classify all optimal codes of length 11 and 12. We finish 
with some concluding remarks in Section IVIII 

II. Codes and Weighted Graphs 

An m-weighted graph is a triple G — (V, E, W), where V 
is a set of vertices, E C V x V is a set of edges, and W 
is a set of weights from GF(m), such that each edge has an 
associated non-zero weight. In an unweighted graph, which 
is simply described by a pair G = (V, E), we can consider 
all edges to have weight one. A graph with n vertices can be 
represented by an n x n adjacency matrix T, where the element 
Tij = W({i,j}) if {i,j} G E, and T^j = otherwise. 
A loop-free undirected graph has a symmetric adjacency 
matrix where all diagonal elements are 0. In a directed graph, 
edges are ordered pairs, and the adjacency matrix is not 
necessarily symmetric. In a colored graph, the set of vertices is 
partitioned into disjoint subsets, where each subset is assigned 
a different color. Two graphs G = (V, E) and G' = (V, E') 
are isomorphic if and only if there exists a permutation 7r 
of V such that {u,v} G E {7r(u), 7r(w)} G E'. 

For weighted graphs, we also require that edge weights are 
preserved, i.e., W({u,v}) = W({tv(u), 7t(u)}). For a colored 
graph, we further require the permutation to preserve the 
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graph coloring, i.e., that all vertices are mapped to vertices 
of the same color. The automorphism group of a graph is 
the set of vertex permutations that map the graph to itself. 
A path is a sequence of vertices, (vx, V2, ■ ■ ■ , Vi), such that 
{i>i, ^2}, {^2, ^3}, . . . , {vi—XjVi} £ E. A graph is connected 
if there is a path from any vertex to any other vertex in the 
graph. 

If an additive code over GF(9) has a generator matrix of 
the form C = T + col, where / is the identity matrix, ui is 
a primitive element of GF(9), and T is the adjacency matrix 
of a loop-free undirected 3-weighted graph, we say that the 
generator matrix is in standard form. A generator matrix in 
standard form must generate a code that is self-dual with 
respect to the Hermitian trace inner product, since it has full 
rank over GF(3) and CC ~ T 2 + T — I only contains entries 
from GF(3), and hence the traces of all elements of lo 2 CC 
will be zero. 

In the context of quantum codes, it was shown by Schlinge- 
mann [18] and by Grassl, Klappenecker, and Rotteler [19] that 
every self-dual additive code is equivalent to a code with a 
generator matrix in standard form. Essentially, the same results 
was also shown by Bouchet l20l in the context of isotropic 
systems. The algorithm given in Fig. Q] can be used to perform 
a mapping from a self-dual additive code to an equivalent code 
in standard form. Note that we can write the generator matrix 
C = A + luB as an n x 2n matrix (A \ B) with elements from 
GF(3). Steps 1 and 2 of the algorithm are used to obtain the 
submatrices A and B. If B now has full rank, we can simply 
perform the basis change B^ 1 {A B) = (T \ I) to obtain 
the standard form. Elements on the diagonal of T can then 
always be set to zero by operations a + bui >—> a + yb + bco, 
for y £ GF(3), corresponding to symplectic matrices (5 j). 
Hence step 12 of the algorithm preserves code equivalence. 
In the case where B has rank k < n, we can assume, after a 
basis change, that the first k rows and columns of B form a 
k x k invertible matrix. This is done in step 5, and the result is 
a permutation of the coordinates of the code. By the operation 
c 1 — y coc, for c = a + bco, corresponding to the symplectic 
we can replace column ai by —bi and bi by 



matrix 
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eti. In this way, we "swap" the n — k last columns of A and B 
in steps 7 and 8. It has been shown that it then follows from 
the self-duality of the code that the new matrix B must have 
full rank |4), ||2TI . and that the matrix T obtained in step 11 
will always be symmetric. 

As an example, consider the (4,3 4 ,3) code generated by 
C which by the described algorithm is transformed into 
the standard form generator matrix C , corresponding to the 
weighted graph depicted in Fig. [2 
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It is known that two self-dual additive codes over GF(4) are 
equivalent if and only if their corresponding graphs are related 
by a sequence of graph operations called local complementa- 
tions (LC) Q, ||20l , II2TI and a permutation of the vertices. We 
have previously used this fact to devise an algorithm to classify 



Require: C generates a self-dual additive code over GF(9). 
Ensure: C generates an equivalent code in standard form. 

U Air- Tr(wC) 
2: B i- Tr(w 2 C) 
3: k <— rank(.B) 
4: if k < n then 

5: Permute rows and columns of B such that the first k 
rows and columns form an invertible matrix. Apply the 
same permutation to the rows and columns of A. 

6: for i = k + 1 to n do 

7: Swap columns and bi 

8: Gtj i di 

9: end for 
10: end if 

ii: r i- b- x a 

12: Set all diagonal elements of T to zero. 

13: C' <- T + LOI 

14: return C 

Fig. 1. Algorithm for Mapping a Code to Standard Form 

1 
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Fig. 2. Graph Representation of the (4, 3 4 ,3) Code 

all self-dual additive codes over GF (4) of length up to 12 [ 1 2 1 . 
The more general result that equivalence classes of self-dual 



additive codes over GF(q 



can be represented as orbits 



of m-weighted graphs with respect to a generalization of LC 
was later shown by Bahramgiri and Beigi [22]. We used this 
to classify all self-dual additive codes over GF(9), GF(16), 
and GF(25) up to lengths 8, 6, and 6, respectively [4|. The 
main obstacle with this approach is that the sizes of the LC 
orbits of weighted graphs quickly get unmanageable as the 
number of vertices increase. We have therefore devised a new 
method for checking code equivalence, which is described in 
the next section. This algorithm uses a graph representation of 
self-dual additive codes over GF(9) that is not related to the 
representation described in this section, and does not require 
the input to be in standard form. However, the weighted graph 
representation will still be very useful for reducing our search 
space. 

III. Equivalence Graphs 

To check whether two self-dual additive codes over GF(9) 
are equivalent, we modify a well-known algorithm used 
for checking equivalence of linear codes, described by 
Ostergard [17|. The idea is to map a code to an unweighted, 
directed, colored equivalence graph such that the automor- 
phism groups of the code and the equivalence graph coincide. 
An important component of the algorithm is to find a suitable 
coordinate graph. For self-dual additive codes over GF(9), we 
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Fig. 3. Coordinate Graphs for Codes over GF(9): (a) Trace-Hermitian Self- 
Dual Additive (b) Linear (c) Hermitian Self-Dual Linear (d) Euclidean Self- 
Dual Linear 



Fig. 4. Example of Equivalence Graph 



need to construct a graph G on eight vertices, labeled with the 
non-zero elements of GF(9), whose automorphism group is 
Sp 2 (3). This graph, shown in Fig. [3] |(a)| was found by adding 
directed edges (<rl, aoS) for all a £ Sp 2 (3). This ensures that 
Sp 2 (3) C Aut(G). We then verified that |Aut(G) = 24 
which implies that Aut(G) = Sp 2 (3). 

Fig. [3] also shows examples of coordinate graphs for some 
other families of codes over GF(9). In the original algorithm 
for checking equivalence of linear codes ifTTl , the coordinate 
graph shown in Fig. [3] |(b)| would be used. This graph has an 
automorphism group of size eight, corresponding to the fact 
that multiplication of a coordinate by any non-zero element 
of GF(9) preserves linearity. The more restrictive coordinate 
graph for Hermitian self-dual linear codes over GF(9) is 
shown in Fig. [3] |(c)| This graph has an automorphism group 
of size four, since only multiplication by x £ GF(9) where 
x 4 = 1 is permitted in this case. Finally, Fig. [3] |(d)| shows 
a graph with automorphism group of size two. This is the 
coordinate graph for Euclidean self-dual linear codes over 
GF(9) where multiplication by ±1 are the only permitted 
operations. Coordinate graphs of this type were used by 
Harada and Ostergard to classify Euclidean self-dual codes 
over GF(5) up to length 16 [23] and over GF(7) up to length 
12 Il24l . 

To construct the equivalence graph of a code, we first add 
n copies of the coordinate graph, each copy representing one 
coordinate of the code. We then need a deterministic way to 
find a set of codewords that generates the code. Taking all 
codewords would suffice, but the following approach yields 



a smaller set and hence a more efficient algorithm. First, 
we check if the set of all codewords of minimum weight 
d generates the code. If it does not, we add all codewords 
of weight d + 1, then all codewords of weight d + 2, etc, 
stopping once we have a set that spans the code. For each 
codeword cj in the resulting set, we add a codeword vertex 
Vi to the equivalence graph. Let the codeword vertices have 
one color, and the other vertices have a different color. Edges 
are added between Vi and the coordinate graphs according 
to the non-zero coordinates of the codeword Ci, e.g., if Cj 
has u in coordinate j, then there is an edge between Vi 
and the vertex labeled to in the jth coordinate graph. As an 
example, Fig. [4] shows the case where c\ = (u)l •1). The 
resulting equivalence graph is finally canonized, i.e., relabeled, 
but with coloring preserved, using the nauty software ll25l . If 
two graphs are isomorphic, their canonical representations are 
guaranteed to be the same. 

Applying a canonical permutation to the vertices of an 
equivalence graph corresponds to permuting the coordinates 
of the corresponding code, applying elements from Sp 2 (3) 
to each coordinate, and sorting the codewords a in some 
canonical order. If two codes are equivalent, their canonical 
equivalence graphs will therefore be identical. Furthermore, 
the automorphism group of a code is equivalent to the auto- 
morphism group of its equivalence graph. This follows from 
the fact that any automorphism of the equivalence graph must 
be one out of 24™n! possibilities, i.e., the n! permutations of 
the n coordinate subgraphs, and the 24 automorphisms from 
Sp 2 (3) of each coordinate subgraph. No other automorphisms 
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are possible. In particular, permuting the codeword vertices 
will never be an automorphism, since all codewords must be 
distinct. Since it is known (2) that coordinate permutations 
and Sp 2 (3) applied to the coordinates of a code preserve its 
weight enumerator, additivity, and self-duality, this must also 
be true for any automorphism of the equivalence graph. 

IV. Classification Algorithm 

We have seen that every weighted graph corresponds to a 
self-dual additive code, and that every self-dual additive code, 
up to equivalence, has a standard form representation as a 
weighted graph. It follows that we only need to consider 3- 
weighted graphs in order to classify all self-dual additive codes 
over GF(9). Permuting the vertices of a graph corresponds 
to permuting coordinates of the associated code, which means 
that we only need to consider these graphs up to isomorphism. 
Moreover, we can restrict our study to connected graphs, 
since a disconnected graph represents a decomposable code. 
A code is decomposable if it can be written as the direct 
sum of two smaller codes. For example, let C be an (n, 3 n , d) 
code and C an (n',3™ , d') code. The direct sum, C © C = 
{u\\v | u £ C,v £ C'}, where || means concatenation, is 
an (n + n', 3™+™ , min{<i, d'}) code. All decomposable codes 
of length n can be generated easily once all indecomposable 
codes of length less than n are known. 

To classify codes of length n, we could take all non- 
isomorphic connected 3-weighted graphs on n vertices, map 
the corresponding codes to equivalence graphs, and canonize 
these. All duplicates would then be removed to obtain one 
representative code from each equivalence class. However, a 
much smaller set of graphs is obtained by taking all possible 
lengthenings ll26ll of all codes of length n — 1. A generator 
matrix in standard form can be lengthened in 3 n_1 — 1 
ways by adding a vertex to the corresponding graph and 
connecting it to all possible combinations of at least one 
of the old vertices, using all possible combinations of edge 
weights. This corresponds to adding a new non-zero row 
r € GF(3) n and column r T to the adjacency matrix, with 
zero in the last coordinate. Only half of the lengthenings 
need to be considered, as adding the row r is equivalent to 
adding r. (Since multiplying the last row and column in the 
corresponding generator matrix by —1 would preserve code 
equivalence.) We have previously shown O, using the theory 
of local complementation of weighted graphs, that the set of 
*n-i — 2~~^ codes obtained by lengthening one representative 
from each of the i n -i equivalence classes of indecomposable 
codes of length n — 1 must contain at least one representative 
from each equivalence class of the indecomposable codes of 
length n. 

Removing possible isomorphisms from the set of length- 
ened graphs, using nauty 11251 . speeds up our classification 
significantly. A set of non-isomorphic graphs that have already 
been processed, as large as memory resources permit, can even 
be stored between iterations, and new graphs can be checked 
for isomorphism against this set. For each graph that is not 
excluded by such an isomorphism check, the corresponding 
code must be mapped to an equivalence graph, as described in 



Require: C n -i contains one graph representation of each 
inequivalent indecomposable code of length n — 1. 

Ensure: C n contains one graph representation of each in- 
equivalent indecomposable code of length n. 

l: C n <- 

2: for all C G C n _i do 

3: E <— 3 " 2 ~ x lengthenings of C 

4: Remove isomorphisms from E 

5: for all E € E do 

6: d<— minimum distance of E 

7: S <— all codewords of weight d from E 

8: while S does not generate E do 

9: d «- d + 1 

10: S <— S U all codewords of weight d from E 

li: end while 

12: Q ■(— equivalence graph given by S 

13: Q' <— canonize Q 

14: G <— graph representation of code given by Q' 

15: if G £ C n then 

16: C n i-C n UG 

17: end if 

18: end for 

19: end for 

20: return C n 

Fig. 5. Classification Algorithm 



Section Hn] The equivalence graph is canonized and compared 
against all previously observed codes, which are stored in 
memory. Since the equivalence graphs will be large, typically 
containing thousands of vertices for n = 10, we map the 
equivalence graph to a canonical generator matrix by taking 
the first n linearly independent codewords corresponding to 
codeword vertices in their canonical ordering. This generator 
matrix can further be mapped to a canonical standard form, 
as described in Section [Til which means that only ( ? 2 l ) ternary 
symbols need to be stored for each code. An outline of the 
steps of our classification algorithm is listed in Fig. [5] 

Note that the special form of a generator matrix in standard 
form makes it easy to find all codewords of low weight, 
which is necessary to construct the equivalence graph. If C 
is generated by C = T + ujI, then any codeword formed by 
taking GF(3)-linear combinations of i rows of C must have 
weight at least i. This means that we can find all codewords of 
weight i by only considering combinations of at most i rows 

of a 



V. Codes of Length 9 and 10 

Using the algorithm described in Section |IV] we have 
classified all self-dual additive codes over GF(9) of length up 
to 10. TableUgives the values of i n , the number of inequivalent 
indecomposable codes of length n, and the values of t n , the 
total number of inequivalent codes of length n. Note that the 
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TABLE I 

Number of Indecomposable (i„) and Total Number (t n ) of 
Self-Dual Additive Codes over GF(9) of Length n 



n 1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


in 1 


1 


1 


3 


5 


21 


73 


659 


17 589 


2 803 404 


tn 1 


2 


3 


7 


13 


39 


121 


817 


18 525 


2 822 779 




numbers t n are easily derived from the numbers i n by using 
the Euler transform [27|: 



Tables [TT] and [TTT] list the numbers of indecomposable codes 
and the total number of codes, respectively, by length 
and minimum distance. In Table [IVj we count the num- 
ber of distinct weight enumerators. There are obviously 
too many codes of length 9 and 10 to list all of them 
here, so an on-line database containing one representative 
from each equivalence class has been made available at 
http://www.ii.uib.no/~larsed/nonbinary/ 

Generator matrices for all the extremal codes of length 9 
and 10 were given in [4]. Our classification confirms that there 
are four extremal (9, 3 9 , 5) codes, all with weight enumerator 
W(y) = 1 + 252y 5 + 1176y 6 + 3672y 7 + 7794y 8 + 678% 9 , 
and with automorphism groups of size 72, 108, 108, and 432, 
and that there is a unique extremal (10, 3 10 , 6) code with 
weight enumerator W(y) = l + 1680y 6 + 2880y 7 + 14040y 8 + 
22160y 9 + 18288y 10 and automorphism group of size 2880. 
The classification of near-extremal codes of length 9 and 10 is 
new. We find that there are 4370 near-extremal (9, 3 9 , 4) codes 
with 25 distinct weight enumerators and 13 different values 
for |Aut(C)|. The weight enumerators that exist are given by 
W 9 , a (y) = l + (4 + 2a) 2 / 4 + (244-4ck)y 5 + (1168-4a) 2 / 6 + 
(3704+16a)2/ 7 +(7766-14a)y 8 +(6796+4a)y 9 for all integer 
values < a < 24. Table [VI gives the number of (9, 3 9 ,4) 
codes for each possible weight enumerator and automorphism 
group size. To highlight a few codes with extreme properties, 
we list generator matrices for the code with automorphism 
group of maximal size (288) and one of the codes with 
weight enumerator Wgfi{y), i.e., with the minimal number of 
minimum weight codewords. For the latter case, we choose the 
unique code with maximum number of automorphisms (16). 
In the following, "-" denotes —1 in generator matrices: 



(uj - 1 1 - - 1 - 1\ 

- oj 1 1 1 1 - 
1 1 uj 1 1 1 - 
1 1 lu 1 1 - 00 

- 1 1 1 u - 1 

- 1 - u 1 
1 1 - OwO 1 

- - 1 1 u 

\ 1 - 1 1 uj J 



L '|Aut|=288 



r^n— 9 

°a=0, |Aut| = 16 



/« 1 1 11 

1 uj 1 - - 1 - 1 

1 1 u 10 1 

1 - - u 1 1 - 1 

1 - - 1 uj 1 1 - 

- 1 - 1 uj - - 1 

- - 1 1 - u 1 - 

- 1 - 1 - 1 u - 
\ - 1 1 - 1 - - uj j 

We find that there are 4577 near-extremal (10, 3 10 , 5) 
codes with 10 distinct weight enumerators and 20 dif- 
ferent values for |Aut(C)|. The weight enumerators that 
exist are given by W Wt0l (y) = 1 + (44 + Aa)y 5 + 
(1460 - 20a)?/ 6 + (3320 + 40a)y 7 + (13600 - 40a)y 8 + 
(22380 + 20a)y 9 + (18244 - 4a)y 10 for integer values a e 
{0, 9, 12, 13, 16, 18, 21, 22, 24, 25}. Table EDgives the number 
of (10, 3 10 , 5) codes for each possible weight enumerator and 
automorphism group size. We give generator matrices for the 
unique codes with automorphism groups of size 2880 and 288, 
as well as the unique code with weight enumerator Wxo,o(y): 



/Willi 
' 1 UJ 1 1 1 



s~in= 10 
L/ |Aut|=2880 



- 

1 - 

1 uj 1 1 1 1 



1 uj 1 - 
1 1 uj - 1 

1 - - OJ 1 



- 1 
1 - 

VlO 



1 

lOllull 
- 1 1 1 uj 1 
1 - 1 1 1 uj 



111 



1 1 uj) 



r^n— 10 
°|Aut|=288 



(uj - 1 1 

Ou - - 1010- 

- uj - 1 1 1 - 
uj 1 - - 1 1 

1 1 1 1^ 
uj 1 - - 1 

1 - lu - 1 
001-0 - - wOl 

1 - 1 1 - 1 uj - 
\ll-10101-w 



(uj 1 - - 

' 1 uj 1 1 

1 oj 1 

- 1 uj 



/~in= 10 



1 1 



0\ 

1 - ' 





- 1 

1 1 1 



- - 1 

1 1 

Vo - o 



uj 



1 - 

1 UJ 1 

- 1 UJ 

- - 1 UJ 

-1-1 

-001 



1 - 

1 

UJ 

1 



That our classification of all codes up to length 10 is correct 
has been verified by the mass formula (fTJ. This required 
us to also calculate the sizes of the automorphism groups 
of all decomposable codes, which was simplified by the 
observation that for a code C = k\C\ © • • • © k m C m , where 
kjCj = ©-ixCj, |Aut(C)| = n™i h\ |Aut(C 4 )| fc \ 

Table I VIII gives the numbers of codes with trivial automor- 
phism group by length and minimum distance. We find that the 
smallest codes with trivial automorphism group are 35 codes 
of length 8. (Note that automorphism group sizes were not cal- 
culated in the previous classification of codes of length 8 (4).) 
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TABLE II 

Number of Indecomposable Self-Dual Additive Codes over GF(9) of Length n and Minimum Distance d 



d\n 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


2 


1 


1 


2 


4 


15 


51 


388 


6240 


418 088 


? 


? 


3 






1 


1 


5 


20 


194 


6975 


893 422 


? 


? 


4 










1 


2 


77 


4370 


1487 316 


? 


7 


5 
















4 


4577 


56005 876 


7 


6 


















1 




6493 


All 


1 


1 


3 


5 


21 


73 


659 


17 589 


2 803 404 


? 


7 



TABLE III 

Total Number of Self-Dual Additive Codes over GF(9) of Length n and Minimum Distance d 



d\n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


1 


1 


1 


2 


3 


7 


13 


39 


121 


817 


18 525 


2 822 779 


? 


2 




1 


1 


3 


5 


20 


60 


424 


6358 


418931 


7 


7 


3 








1 


1 


5 


20 


195 


6976 


893 429 


7 


7 


4 












1 


2 


77 


4370 


1487 316 


7 


7 


5 


















4 


4577 


56 005 876 


7 


6 




















1 




6493 


All 


1 


2 


3 


7 


13 


39 


121 


817 


18 525 


2 822 779 


> 2 30 


> 2 41 



TABLE IV 

Number of Distinct Weight Enumerators of Indecomposable Codes of Length n and Minimum Distance d 



d\n 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


2 


1 


1 


2 


4 


14 


42 


202 


1021 


8396 


7 


? 


3 






1 


1 


3 


9 


33 


170 


1133 


7 


? 


4 










1 


1 


9 


25 


345 


7 


? 


5 
















1 


10 


48 


7 


6 


















1 




27 


All 


1 


1 


3 


5 


18 


52 


244 


1217 


9885 


7 


7 



We give the generator matrix for one (8, 3 8 , 4) code with trivial 
automorphism group. Generator matrices for the other codes 
can be obtained from http://www.ii.uib.no/~larsed/nonbinary/ 



/u 
oj - 
- LJ 



/^m—S 
°|Aut|=2 



V 



11 

111- 

M-10-0 

- w - 1 
1 u 1 
0-1^00 

- oj - 
0100-u 



We observe that codes with minimum distance d < 2 always 
have nontrivial automorphisms, and this can be proved as 
follows. For d = 1, we can assume that the first row of 
a standard form generator matrix is (wO- ■•()). Then (j;?) 
applied to the first coordinate of the code is an automorphism 
of order 3. Multiplying the first coordinate by —1 has the same 
effect as multiplying the first row by —1 and is therefore an 
automorphism of order 2. Including the trivial automorphism, 
we have that |Aut| > 12. There are codes of length 9 with 
d = 1 and |Aut| = 12 which shows that the bound is tight. 
For d = 2, we can assume that the first row of a standard 
form generator matrix is (wlO • ■ - 0). Then (J §) applied to 



the first coordinate and ( J i ) applied to the second coordinate 
of the code has the same effect as adding the first row of 
the generator matrix to the second row, and is hence an 
automorphism of order 3. Swapping the first two coordinates 
is an automorphism of order 2, since it has the same effect as 
the following procedure: Add the first row to itself, then add 
the first row to each row i > 2 where the value in position i 
of the second column is 1 , and add twice the first row to each 
row i > 2 where the value in position i of the second column 
is 2. Finally apply ( ° § ) to the first column, and to the 

second column. Again, including the trivial automorphism we 
get the bound |Aut| > 12, and the existence of codes of length 
8 with d = 2 and Aut =12 proves that the bound is tight. 

VI. Optimal Codes of Length 1 1 and 12 

When we lengthen an (n,3 n ,d) code, as described in 
Section |IV] we always obtain an (n+ 1, 3 n+1 , d') code where 
d' < d + 1 ll26l . It follows that given a classification of all 
codes of length n and minimum distance at least d, we can 
classify all codes of length n+l and minimum distance at least 
d+1. There are no (11, 3 , 6) codes, but by lengthening the 
1491 894 (10,3 10 ,d) codes for d > 4, we are able to obtain 
all optimal (11, 3 11 , 5) codes. To quickly exclude codes with 



TABLE V 

Number of (9, 3 9 , 4) Codes with Weight Enumerator Wg iCX (y) and |Aut(C)| = /3 





2 


4 


6 


8 


12 


16 


24 


32 


36 


48 


72 


144 


288 


All 







3 




1 




1 
















5 


1 


2 


1 


2 






















5 


2 


15 


21 




4 




















40 


3 


15 


13 


1 


3 


2 




1 














35 


4 


125 


52 






12 


2 


2 






2 








195 


5 
6 


85 

338 


8 
93 




11 


2 


1 
















93 
445 


7 
8 


165 
561 


53 
150 


2 


9 
11 


2 


2 


2 


1 






1 






236 
723 


9 


173 


20 


6 


7 






1 














207 


10 


522 


154 


4 


7 


7 




2 






2 








698 


11 


157 


53 




15 




















225 


12 


356 


143 


2 


4 


3 


2 
















510 


13 


119 


25 


2 


6 




















152 


14 


229 


114 




11 




2 
















356 


15 


42 


28 


1 


16 


1 


2 
















90 


16 


96 


62 




8 


2 


3 


2 




4 






2 


1 


180 


1 / 


1 c 

ID 


n 

y 




t: 





















JO 


18 


23 


33 




2 


1 


6 




2 












67 


19 


9 


4 




6 




2 


2 














23 


20 


8 


23 




6 








2 












39 


21 




2 


2 








1 














5 


22 


1 


3 








2 




1 












7 


23 








1 




















1 


24 
















3 












3 


All 


3056 


1067 


22 


134 


32 


25 


13 


9 


4 


4 


1 


2 


1 


4370 



TABLE VI 

Number of (10, 3 10 , 5) Codes with Weight Enumerator Ww, a (y) and | Aut(C) | = /3 



a\f3 


2 


4 


6 


8 


10 12 


16 


20 


24 


32 


36 


40 48 


64 72 144 


192 240 288 


2880 All 





























1 


1 


9 












1 










1 




1 1 


4 


12 












2 
















2 


13 




3 








3 






1 










7 


16 


10 


5 




2 












1 








18 


18 


30 


24 


4 


4 


8 


4 
















74 


21 


190 


77 


2 


20 


2 


2 




4 








3 




300 


22 


467 


72 




4 




1 
















544 


24 


2321 


172 


4 


4 


1 5 


1 










1 1 






2510 


25 


777 


247 


12 


39 


14 


10 


3 


2 


1 


2 


2 


4 1 1 


1 


1 1117 



All 3795 600 22 73 1 29 24 3 6 2 3 4 1 7 1 1 1 2 1 1 4577 
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TABLE VII 

Number of Codes of Length n and Minimum Distance d with 
Trivial Automorphism Group 



d\n 


< 7 


8 


9 


10 


11 


12 


< 2 




















3 





32 


4518 


832 878 


7 


7 


4 





3 


3056 


1419 861 


? 


7 


5 









3795 


55 865 753 


7 


6 













3445 


All 





35 


7574 


2 256 534 


7 


7 



d < 5, we checked the minimum distance of each lengthened 
code before checking for code equivalence in this search. 

We find that there are 56005 876 optimal (11, 3 11 , 5) codes 
with 48 distinct weight enumerators and 24 different values 
for |Aut(C)|. The weight enumerators that exist are given by 
W Ut a(y) = 1 + (12 + 2a)y 5 + (888 - 6a)y 6 + 3960?/ + 
(14970 + 20a)y 8 + (42500 - 30a)y 9 + (66240 + 18a)y 10 + 
(48756 — 4a)?; 11 for all integer values 6 < a < 50 as well 
as a £ {0, 54, 60}. Observe that the number of codewords of 
weight 7 is constant for all codes. Table IVIHl gives the number 
of (11, 3 11 , 5) codes for each possible weight enumerator and 
automorphism group size. We give generator matrices for the 
unique codes with automorphism group of size 47520 and 
1440, as well as the unique code with weight enumerator 
Wu,o(y): 



We find that there are 6493 optimal (12, 3 12 , 6) codes with 
27 distinct weight enumerators and 32 different values for 
Aut(C)|. The weight enumerators that exist are given by 
Wi2, a (y) = 1 + (480 + 4a)y 6 + (3456 - 2Aa)y 7 + (15120 + 
60ajy 8 + (55520 - 80a)y 9 + (133920 + 60a)y 10 + (19536 - 
24a)?/ 11 + (129408 + 4a)y 12 for all integer values a G {0, 
1, 3, 4, 7, 9, 12, 13, 16, 19, 21, 25, 27, 28, 31, 36, 37, 
39, 43, 48, 49, 52, 57, 63, 64, 81, 144}. Table |IX] gives 
the number of (12, 3 12 , 6) codes for each possible weight 
enumerator and automorphism group size. Generator matrices 
for all the optimal codes of length 12 can be obtained from 
http://www.ii.uib.no/~larsed/nonbinary/ We here list generator 
matrices for the unique code with maximal automorphism 
group size (2 280960) and a code with weight enumerator 
Wi2 t o (y). In the latter case, we choose the single code with 
maximal number of automorphisms (11 520). 



^m=12 

°|Aut|=2 280 960 



1 1 - - 0\ 

1 - - 1 

- - 1 1 
1-1-0 

- 1 - 1 

uj 

1 - - Owl 1 1 1 - 

- - 1 1 1 uj 1 1 1 - 

- 1 - 1 1 w 1 1 - 

oil-omul- 

10-lOllllu- 
\0 uj ) 



UJ - - - - 

- UJ - - - 

- - UJ - - 

- - - UJ - 

- - - - UJ 



L/ |Aut|=47 520 



c\ 



Aut| = 1440 



ill - 11 - 11 

1 uj 1 1 

- 1 uj - - 1 
lOOwOll- 
1 1 



- - 
1 - - 
1 - 1 
111 



0^0001 
1 uj 1 
100w0- 

- uj 

- 1 1 - uj 
- - - 



1 0\ 
10 1* 

111 

- 1 

- 

- 


- 



10 
\ 1 1 1 



1 - 1 - 



(UJ 

1«-11- 

- - ulOl 
1 1 1 uj - 

- lOOul 
- 1 - 1 uj 

a; 

QQ--110w001 

01- 01-00^0- 
1 - - 1 uj - 

\1 - 1 1 - 1 - - w 




--011 

- 

- 

1 - 

- 1 



r~in— 11 



(UJ 

1 


1 




1 - 

uj I 

1 uj 1 

1 u 

1 - 





1 



- 1 0\ 
1 



- 1 

10 0-0 



1-0 

\ - 1 




- - 
-000 

1 





1 



- uj - 

- uj - - 

UJ I 

1 - 
1 - 



1 uj Q 
10 uj) 



(uj 

- uj - - - 

- - uj - - 

- - - UJ - 

- - - - UJ 



(~<n=12 

°a=0, |Aut|=ll 520 



1 - 

- 



1 

1 - 



- 1 0\ 
1-1 

-11- 

-0-1 
1- 



uOOOOOO 

IOOIIOcj 11 



1 





0-0 
--00 
10 
1-10 
-1-0 



uJ 



1 

- - UJ - I 

UJ 1 

1111a; 
1111 



1 
1 
1 
1 

1 uj) 



VII. Conclusion 

According to the mass formula bound ©, the total number 
of codes of length 11 and 12 are t n > 1592 385 579 and 
ti2 > 2 938 404 780 748, which makes complete classifica- 
tions infeasible, at least with our computational resources. 
Running our algorithm on a typical desktop computer, the 
classification of codes of length n was completed in less 
than five minutes for n < 8, about two hours for n = 9, 
and about a week for n — 10. Most of this time is spent 
canonizing the equivalence graphs with nauty, and far more 
time is used on codes with large automorphism groups than 
on codes with trivial or small automorphism groups. This 
means that our previous classification algorithm [4|, using 
local complementation, might still be useful in some cases, 
since we observe that graphs corresponding to codes with 
large automorphism groups typically have small LC orbits. 
For instance, we could speed up our classification algorithm 
by not only removing isomorphisms from the set of lengthened 
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TABLE VIII 









Number of (11, 3 


11 , 5) Codes wr 


rH Weig 


ht Enumerator Wh. 


AND 


|Aut(C)| = P 




a\p 


2 


4 


6 


8 


10 12 


16 


18 


20 24 


32 36 


40 44 48 72 108 


120 


144 


288 360 432 1440 47 520 


All 



6 






4 








1 
1 




3 


2 








1 
10 


7 
8 


10 


4 
5 




1 








2 












7 
15 


9 


36 


22 


4 


4 


4 








2 










72 


10 


35 


16 


2 


2 


2 


















57 


11 


286 


62 




2 




















350 


12 


217 


37 


15 


7 


2 


1 
















279 


13 


1515 


170 


6 


8 


4 






2 












1705 


14 


1140 


139 




4 




















1283 


15 


7412 


414 


10 


20 


14 






1 


8 


5 








7884 


16 


5234 


192 


4 


13 




2 
















5445 


17 


30 825 


906 




28 




















31759 


18 


19468 


623 


17 


14 


12 


















20134 


19 


108 109 


1606 


26 


24 


8 


















109773 


20 


62 364 


641 




7 


5 


















63 017 


21 


314156 


2701 


16 


42 


27 






2 












316944 


22 


169270 


1928 


4 


30 


10 


















171242 


23 


780271 


4123 


5 


40 




















784439 


24 


385400 


1508 


38 


32 


1 8 


2 




1 


2 3 


1 




1 


1 1 


386999 


25 


1649942 


5666 


42 


33 


2 


















1 655 685 


26 


754 931 


4249 




44 


10 


3 




4 












759 241 


27 


2990527 


7882 


61 


36 


44 




2 




2 


6 








2 998560 


28 


1 266 193 


2610 


20 


19 


6 








1 


2 








1268 851 


29 


4671482 


9256 


18 


36 












2 








4 680794 


30 


1832 724 


6641 


41 


50 


20 


















1 839476 


31 


6241 827 


10336 


98 


39 


10 


















6 252 310 


32 


2266449 


3048 




45 




6 
















2 269 548 


33 


7110043 


10986 


89 


27 


47 








7 


4 2 








7 121 205 


34 


2 377 017 


7970 


44 


66 


4 


6 
















2 385 107 


35 


6821413 


10684 


6 


22 




















6 832125 


36 


2084454 


3159 


46 


29 


8 


















2087 696 


37 


5 388 851 


9356 


99 


22 


6 


















5 398334 


38 


1 475 547 


6545 




30 


6 


















1482128 


39 


3403 383 


7317 


65 


27 


2 36 












1 






3 410 831 


40 


810399 


2084 


34 


48 


2 


4 








2 








812 573 


41 


1 645 374 


5231 




13 




















1650618 


42 


334 536 


3308 


35 


39 


28 








3 


1 2 




1 




337 953 


43 


579 338 


2764 


32 


6 


6 






2 












582 148 


44 


94 833 


664 




10 




2 
















95 509 


45 


137 174 


1487 


18 




26 








2 










138 707 


46 


21 818 


713 


4 




4 


















22 539 


47 


18 178 


353 
























18531 


48 


1901 


113 


12 


2 


2 


2 




1 












2033 


49 


1275 


174 


3 




1 4 






2 












1459 


50 


392 


80 




4 




















476 


54 


4 


9 




4 




















17 


60 










2 










2 






1 1 1 


7 



All 55865753 137782 918 929 4 369 28 4 4 13 3 30 4 6 4 13 4 1 2 1 1 1 1 1 56005876 
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TABLE IX 

Number of (12, 3 12 , 6) Codes with [Aut(C)| = p and Weight Enumerator Wi2, a (y) 







1 


a 

j 


4 


7 


Q 


1 2 


1 ^ 
i j 


i fx 

1 u 


1 

i j 


9 | 


9S 

Z J 


97 


~>8 
— O 


J 1 


3fi 


3 / 


tO 

3 y 


zR 
4 J 


-fo 




S9 

3Z 


57 




U4 


8 1 


144 


All 


9 

Z 


J J 


1 1 7 


34 


IzU 


1 8A 
1 oD 


7 MO 

zuy 


1 ^8 

1 Jo 


33 O 


T9^ 

3Z3 


a -i 8 

44 o 


41 o 


Z30 


1 AH 
1DU 


9A8 
ZOo 


1 OZ 


^7 


8A 
oO 


^a 

30 


O 




O 














1AA ^ 
3443 


A 
4 


OZ 


1 j J 


1 uz 


1 zT7 
14 / 


1 HA 
1 o4 


8S 


1 Z4 


9 1 A 
Z 1 4 


1 O 1 


1 88 
1 oo 


999 

zzz 


1 1 3 


on 


1 34 


88 
00 


3D 


4o 


30 


9zl 
Z4 


7 

1 


1 O 














99^ 
ZZ30 







c 

J 


z 


o 


Z 


o 






A 
4 




g 
o 


9 
Z 


zl 

4 




9 

z 


4 










1 
1 














40 


Q 




Zj 


Zj 


Q 


2 1 

3 1 


4U 


10 


JZ 


Zo 


Zo 


ZZ 


1(5 


1 -1 
14 


1 -1 
14 


ZD 


10 


1 3 





A 
4 






g 
O 


9 

z 












3 /0 


1 Z 


J 





1 Z 


1 u 


4 


iy 


1 z 


Q 
O 


"7 
/ 





1 Z 


1 





Z 










9 

z 


9 

z 




-> 

Z 


z 












1 ZO 


1 O 


90 
ZU 




Z 


zl 

4 


D 


z 


J 


A 
4 


1 o 


4 


4 


9 
Z 


9 

z 








-1 
J 


9 

Z 


9 
z 




A 
4 
















09 

VZ 


OA 
Z4 


z 


A 
4 


7 

1 


o 




g 
o 






1 

1 




g 
o 




9 

Z 


4 




A 
4 




V 




9 

Z 


9 

Z 






zl 
4 










^zl 
D4 


3Z 


J 






J 






1 
1 
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codes, but also generating and storing a limited number of 
LC orbit members of each graph, and checking new graphs 
for isomorphism against this set. Finding all optimal codes 
of length 11 and 12 required 80 and 320 days of CPU time, 
respectively, and a parallel cluster computer was used for this 
search. We observed that most of this time was spent on 
computing minimum distance to eliminate non-optimal codes, 
and much less time on canonizing the optimal codes. 

Although this paper has focused on codes over GF(9), 
our classification algorithm can be generalized to Hermitian 
self-dual additive codes over GF(q — to 2 ) for any prime 
power m. (One simply needs to find an appropriate coordinate 
graph, as discussed in Section [Til]) The results in this paper 
also has applications beyond the study of additive codes. The 
correspondence between self-dual additive codes over GF(9) 
and 3-weighted graphs means that we have also classified 



particular classes of 3-weighted graphs that should be of 
interest in graph theory. An equivalence class of self-dual 
additive codes over GF(9) maps to an orbit of graphs under 
generalized local complementation H, 1221 . Orbits of graphs 
with respect to local complementation has a long history 
in combinatorics [20|, |28|, with several applications, for 
instance in the theory of interlace polynomials ll29ll . [30]. 
The generalization to weighted graphs is a natural next step. 
The results in this paper also have applications in the field 
of quantum information theory. Our previous classification of 
codes over GF(4) ifPIl has since led to new results in the 
study of the entanglement of quantum graph states [31], and 
the new data obtained in this paper will yield similar insights 
into the properties of ternary quantum graph states. 
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